LOCAL MULTIFRACTAL ANALYSIS 



JULIEN BARRAL, ARNAUD DURAND, STEPHANE JAFFARD AND STEPHANE SEURET 

Abstract. We introduce a local multifractal formalism adapted to functions, 
measures or distributions which display multifractal characteristics that can 
change with time, or location. We develop this formalism in a general frame- 
work and we work out several examples of measures and functions where this 
setting is relevant. 



1. Introduction 

Let / denote a function, a positive Radon measure, or, more generally, a distri- 
bution defined on a nonempty open set fi. One often associates with / a pointwise 
exponent, denoted by hf(x), which allows to quantify the local smoothness of / at 
x. On the mathematical side, the purpose of multifractal analysis is to determine 
the fractal dimension of the level sets of the function x h- > hf(x). Let 

E H = {x: h f {x) = H}. 

The multifractal spectrum of / (associated with the regularity exponent ht) is 

d f (H)= dim E„ 

(where dim denotes the Hausdorff dimension, see Defirrition[3j). Multifractal spectra 
yield a description of the local singularities of the function, or measure, under 
consideration. 

Regularity exponents (and therefore the multifractal spectrum) of many func- 
tions, stochastic processes, or measures used in modeling can be theoretically de- 
termined directly from the definition. However, usually, one cannot recover these 
results numerically on simulations, because the exponents thus obtained turn out 
to be extremely erratic, everywhere discontinuous functions it is for instance the 
case of Levy processes [32J, or of multiplicative cascades (see the book [TU , and, 
in particular the review paper by J. Barral, A. Fan, and J. Peyriere) so that a 
direct determination of hf(x) leads to totally instable computations. A fortiori, 
the estimation of the multifractal spectrum from its definition is unfeasable. The 
multifractal formalism is a tentative way to bypass the intermediate step of the 
determination of the pointwise exponent, by relating the multifractal spectrum 
directly with averaged quantities that are effectively computable on experimental 
data. Such quantities can usually be interpreted as global regularity indices. For 
instance, the first one historically used in the function setting (C/(p), referred to as 
Kolmogorov scaling function) can be defined as follows; for the sake of simplicity, 
we only consider in this introduction the function setting and we assume here that 
the functions considered are defined on the whole M. d . 
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Recall that the Lipschitz spaces are defined, for s <G (0,1), and p £ [l,oo], by 
/ G Lip(s, LP(R d )) if / G LP and 3C > such that V7i > 0, 

||/(a: + /i)-/(aO|| L p<C7i'. (1) 

(the definition for larger s requires the use of higher order differences, and the 
extension to p < 1 requires to replace Lebesgue spaces by Hardy spaces, see [51]). 
Then 

C/(p) = P ■ sup{ S : / G Lip( S , L p (R d ))}. (2) 

Initially introduced by U. Frisch and G. Parisi in the mid 80s, the purpose of 
multifractal analysis is to investigate the relationships between the pointwise regu- 
larity information supplied by df(H) and the global regularity information supplied 
by C/(p)- Note that these quantities can be computed on the whole domain of defi- 
nition fi of /, or can be restricted to an open sub-domain uj C f2. A natural question 
is to understand how they depend on the region uj where they are computed. It is 
remarkable that, in many situations, there is no dependency at all on u>; we will 
then say that the corresponding quantity is homogeneous. It is the case for several 
classes of stochastic processes. For instance, sample paths of Levy processes (and 
fields) [24, 25, 33], Levy processes in multifractal time [15], and fractional Brownian 
motions (FBM) almost surely have homogeneous Holder spectra, and, in the case 
of FBM, the Legendre spectrum also is homogeneous, see [33J 123 . In the random 
setting, it is also the case for many examples of multiplicative cascades, see [T2"] . 
Many deterministic functions or measures also are homogeneous (homogeneity is 
usually not explicitly stated as such in the corresponding papers, but is implicit in 
the determination of the spectra). This is for instance the case for self-similar or 
self-conformal measures when one assumes the so-called open set condition, or for 
Gibbs measures on conformal repellers (see for instance [UJ 051 [50] ) • It is also the 
case for many applications, for instance the Legendre spectra raising from natural 
experiments (such as turbulence, see [2] 0] and references therein) are found to be 
homogeneous. 

On the opposite, many natural objects, either theoretical or coming from real 
data, have been shown to be non-homogeneous : Their multifractal characteristics 
depends on the domain f2 over which they are observed: 



It is the case of some classes of Markov processes, see [5] and Section 5.2 
and also of some Markov cascades studied in 

Some self-similar measures when the open set condition is relaxed into the 
weak-separation condition may satisfy the multifractal formalism only when 
restricted to some intervals (see [501 1301 15T1 12"0] ) . 

In applications, many types of signals, which have a human origin, can 
have multifractal characteristics that change with time: A typical exam- 
ple is supplied by finance data, see [2], where changes can be attributed 
to outside phenomena such as political events, but also to the increasing 
sophistication of financial tools, which may lead to instabilities (financial 
crises) and implies that some characteristic features of the data, possibly 
captured by multifractal analysis, evolve with time. This situation is also 
natural in image analysis because of the occlusion phenomenon; indeed, a 
natural image is a patchwork of textures with different characteristics, so 
that its global spectrum of singularities reflects the multifractal nature of 
each component, and also of the boundaries (which may also be fractal) 
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where discontinuities appear. Note that the notion of local Hausdorff di- 
mension which plays a central role in this section, has been introduced in 
|38j precisely with the motivation of image analysis. 
• Functions spaces with varying smoothness have been introduced motivated 
by the study of the relationship between general pseudo-differential oper- 
ators and later by questions arising in PDEs, see |52j for a review on the 
subject; scaling functions with characteristics depending on the location 
are then the natural tool to measure optimal regularity in this context. We 
will investigate this relationship in Section [7] 

This paper will provide new examples of multifractal characteristics which de- 
pend on the domain of observation. In such situations, the determination of a local 
spectrum of singularities for each "component" w C will carry more information 
than the knowledge of the "global" one only. A natural question is to understand 
how the different quantities which we have introduced depend on the region uj where 
they are computed. 

Some of the notions studied in this paper have been already introduced in [S]; 
let us also mention that a local i 9 -spectrum was already introduced in [39] , where 
the authors studied this notion for measures in doubling metric spaces (as well 
as the notion of local homogeneity) and obtained, for instance, upper bounds for 
the dimensions of the sets of points with given lower and upper local dimensions 
using this local concepts. The goal of their approach was to investigate conical 
density and porosity questions. In our paper, on top of measures, we also deal 
with functions, get comparable upper bounds for the multifractal spectra, and the 
examples we develop are very different. 

Let us now make precise the notion we started with, namely pointwise regularity. 
The two most widely used exponents are the pointwise Holder exponent of functions 
and the local dimension of measures. In the following, B(xo,r) denotes the open 
ball of center xq and radius r. 

Definition 1. Let [ibe a positive Radon measure defined on an open subset C K d . 
Let xq € and let a > 0. The measure \i belongs to h a (xo) if 

3C, R > 0, Vr < R, fi(B{x ,r)) < Cr a . (3) 

Let Xq belong to the support of /j,. The lower local dimension of fi at xq is 

j / \ r 7 a/ w }■ r log fj,(B(x ,r)) 

hfj,{xo) — supja : \i € h (xo)} = hmmf . (4) 

r— >0+ log V 

We now turn to the case of locally bounded functions. In this setting, the notion 
corresponding to the lower local dimension is the pointwise Holder regularity. 

Definition 2. Let xq g M. d and let a > 0. Let f : 57 — > M be a locally bounded 
function; f belongs to C a (xo) if there exist C, R > and a polynomial P of degree 
at most a such that 

if \x-x \<R, then \f(x)-P{x-x Q )\<C\x-x \ a . (5) 

The Holder exponent of f at xq is 

ft / M=sup{a:/eC a M}. (6) 



This paper is organized as follows: 
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In Section [2] we recall the notions of dimensions that we will use (both in the 
global and local case), we prove some basic results concerning the notion of local 
Hausdorff dimension, and we recall the the wavelet characterization of pointwise 
Holder regularity. 

In Section|3]we recall the multifractal formalism on a domain in a general abstract 
form which is adapted both to the function and the measure setting; then the 
corresponding version of local multifractal formalism is obtained, and we draw its 
relationship with the notion of germ space. 

In Section [4j we investigate more precisely the local multifractal analysis of 
measures, providing natural and new examples where this notion indeed contains 
more information than the single multifractal spectrum. In particular, we introduce 
new cascade models the local characteristics of which change smoothly with the 
location; here again, we show that the local tools introduced in Section [2] yield the 
exact multifractal characteristics of these cascades. 

In Section [5j we review the results concerning some Markov processes which do 
not have stationary increments; then we show that the notion of local spectrum 
allows to recover the exact pointwise behavior of the Multifractional Brownian 
Motion (in contradistinction with the usual "global" multifractal formalism). 

In Section [6] we consider other regularity exponents characterized by dyadic 
families, and show how they can be characterized in a similar way as the previous 
ones, by log-log plot regressions of quantities defined on the dyadic cubes. 

Finally, in Section [7] the relationship between the local scaling function and 
function spaces with varying smoothness is developed. 

2. Properties of the local Hausdorff dimension and the local 

multifractal spectrum 

2.1. Some notations and recalls. In order to make precise the different notions 
of multifractal spectra, we need to recall the notion of dimension which will be 
used. 

Definition 3. Let AdM. d . If e > and 5 £ [0, d], we denote 

Mf = inf (El^l 5 

where R is an e-covering of A, i.e. a covering of A by bounded sets of 
diameters \Ai\ < e. The infimum is therefore taken on all e-coverings R. 
For any S £ [0, d], the 6 -dimensional Hausdorff measure of A is 

mess (A) = lim Mf . 

There exists 6q £ [0, d] such that 

V < S < So, mesg(A) = +oo and V<5 > <5o, mesg(A) = 0; 

this critical value Sq is called the Hausdorff dimension of A, and is denoted by 
dim(A). By convention, we set dim(0) = — oo. 

In practice, obtaining lower bounds for the Hausdorff dimension directly from the 
definition involve considering all possible coverings of the set, and is therefore not 
practical. One rather uses the mass distribution principle which involves instead 
the construction of a well-adapted measure. 
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Proposition 1. Let A cR d and let /i be a Radon measure such that fi(A) > 0; 

if Vx £ A, limsup /i(5(x,r)) < C then H S (A) > 

r->o r s c 

We will see in Section [3~2l a local version of this result. 

Apart from the Hausdorff dimension, we will also need another notion of dimen- 
sion: The packing dimension which was introduced by C. Tricot, see |58j : 

Definition 4. Let A be a bounded subset o/R d ; if e > 0, we denote by N e (A) the 
smallest number of sets of radius e required to cover A. The lower box dimension 
of A is 

dim R M) = liminf — - , e ^ — -. 

e^o -loge 

The packing dimension of a set A C K d is 

dim p (A) = inf \ sup dim B A: : A C I J Ai ) > (7) 

l i6N V i=i J) 

(the infimum is taken over all possible partitions of A into a countable collection 
AJ. 

2.2. Local Hausdorff dimension. 

In situations where the spectra are not homogeneous, the purpose of multifractal 
analysis is to understand how they change with the location where they are consid- 
ered. In the case of the multiractal spectrum, this amounts to determine how the 
Hausdorff dimension of the set Ef(H) changes locally. This can be performed using 
the notion of local Hausdorff dimension, which can be traced back to |38j (see also 
[8] where this notion is shown to be fitted to the study of deranged Cantor sets). 

Definition 5. Let A C K d ; and x £ K d . The local Hausdorff dimension of A at x 
is the function defined by 

Vx G A, 6im(A,x) — lim dimMn B(x, r)). (8) 

Remarks: 

• The limit exists because, if fli C O2 , then dim(Oi) < dim(fl2); therefore 
the right-hand side of ([8]), being a non-negative increasing function of r, 
has a limit when r — > 0. 

• We can also conisder this quantity as defined on the whole K rf , in which 
case, it takes the value —00 outside of A. 

• The same definition allows to define a local dimension, associated with any 
other definition of fractional dimension; one gets for instance a notion of 
local packing dimension. 

The following result shows that the local Hausdorff dimension encapsulates all 
the information concerning the Hausdorff dimensions of the sets of the form A n uj, 
for any open set to. 

Proposition 2. Let A C M d ; then for any open set lu which intersects A, 

dim(^4 n oS) = sup dim(A, a;). (9) 
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Proof. For r small enough, B r c w, it follows that 

Vie G w, dim(A, x) < dim(A n cj), 

and therefore sup dim( A, x) < dim(A n w). 

Let us now prove the converse inequality. Let (K n ) ne jq be an increasing sequence 
of compact sets such that UK n = cj; then 

dim(Ano;)= lim dim(A n K„). 

71— >00 

Let 5 > be given; then 

Va; G K n , 3r(x) > 0, dim(A n B{x, r)) - dim(^4, x) < S. 

We extract a finite covering of K n from the collection {B(x, r(x))} xeK which yields 
a finite number of points x\, ■ ■ ■ x^ G uj such that K n C (J B(xi, r(xi)); thus 

dim(j4ni<r„) < sup dim(A n B(xj, r(xj))) 

i=l,- ,JV 

< sup dim(A, Xi) + S < sup dim(A, x) + 5. 

Taking S — > and N — > oo yields the required estimate. □ 

Proposition [2] implies the following regularity for the local Hausdorff dimension. 

Corollary 1. Let A be a given subset ofM. d ; then the function x — > dim(^4,a;) is 
upper semi-continuous. 

Proof. We have 

dim(j4,x) = lim dim(A D B(x, r)) = lim sup dim(A, y) =limsup dim(^4,y). 

r->0 r^O yeB (x,r) ' y->x 

□ 

2.3. Wavelets and wavelet leaders. In Section [3] we will describe a general 
framework for deriving a multifractal formalism adapted to pointwise regularity 
exponents. The key property of these exponents that we will need is that they are 
derived from log-log plot regressions of quantities defined on the dyadic cubes. Let 
us first check that it is the case for the pointwise exponent of measures. 
Recall that a dyadic cube of scale j € Z is of the form 



A = 



fci fei + 1 
2T' V 



kd k d + l 
23 ' 1? 



(10) 



where k = (fc 1; . . . kj) G 7L d . Each point Xq G K d is contained in a unique dyadic 
cube of scale j, denoted by Xj(xo). 

Let 3Aj(xo) denote the cube with the same center as Xj(xo) and three times 
wider; it is easy to check that ^ and Q can be rewritten as 

M 0J j^+oo log2"J 

We now show that the Holder exponent of a function can be recovered in a similar 
way, from quantities derived from wavelet coefficients. Recall that orthonormal 
wavelet bases on M. d are of the following form: There exist a function ip and 2 d — 
1 functions V'^ with the following properties: The <p(x — k) (k G Z d ) and the 
2«/2^W(2J'a; - k) (k G Z d , j G Z) form an orthonormal basis of L 2 (R d ). This basis 



LOCAL MULTIFRACTAL ANALYSIS 7 

is r-smooth if ip and the tp^ are C r and if the d a ip, and the d a (pip^\ for \a\ < r, 
have fast decay. Therefore, V/ € L 2 , 

oo 

f( X )= j2^(x-k)+j2J2J2 c l^ {l) ^x~ky 1 (ii) 

the Cfe and c*- fc are the wavelet coefficients of /: 



2*/ /(x)^ W (2 J x- fc)dar, and c k = / f{x)cp(x - k)dx. (12) 



Note that (12) makes sense even if / does not belong to L ; indeed, when using 
smooth enough wavelets, these formulas can be interpreted as a duality product 
between smooth functions (the wavelets) and distributions. 

Instead of the three indices k), wavelets will be indexed by dyadic cubes as 
follows: Since the wavelet index i takes 2 d — 1 values, we can assume that it takes 
values in {0, l} d - (0,...,0); we will use the notations 

d 



2 J 23 + 1 



0,-^) > c A =c}, fc! ^ x {x) = ^{Vx-k). 



Note that the cube A which indexes the wavelet gives information about its loca- 
tion and scale; if one uses compactly supported wavelets, then 3C > such that 
supp(^a) C C ■ A. 

Finally, Aj will denote the set of dyadic cubes A which index a wavelet of scale 
j, i.e. wavelets of the form ?p\(x) — ^ % ^{2?x — k) (note that Aj is a subset of the 
dyadic cubes of side 2 3+1 ). We take for norm on ~R d 

if x = (xi,...,x d ), \x\= sup \xi\; 

i=l,...,d 

so that the diameter of a dyadic cube of side 2~ 3 is exactly 2~' J . 

In the following, when dealing with Holder regularity of functions, we will always 
assume that, if a function / is defined on an unbounded set SI, then it has slow 
increase, i.e. it satisfies 

3C,N>0 \f(x)\<C(l + \x\) N ; 

and, if O R d , then the wavelet basis used is compactly supported, so that, if 
xq € Q, then the wavelet coefficients "close" to xo are well defined for j large 
enough. 

Let / be a locally bounded function, with slow increase. The pointwise Holder 
regularity of / is characterized in terms of the wavelet leaders of /: 

d x = sup \cy\. (13) 

A'C3A 

The assumptions we made on / imply that wavelet leaders are well defined and 
finite. 

We note dj(xo) — d\.t Xo y The following result allows to characterize the Holder 
exponent by the decay rate of the cJa-(x ) when j — > +oo, see [33]. 

Proposition 3. Let a > and let ip\ be an orthonormal basis with regularity 
r > a. If there exists e > such that f £ C e (Q), then 

Vx , fe / (xo) = liminf 1 ° grfA ff , (14) 

j— S.+OG log Z J 
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Hence, the pointwise Holder exponent can be computed from a dyadic family. 
This is also the case for the lower dimension of a measure /i. Indeed, if 3\j(xo) 
stands for the dyadic cube A of generation j neighboring Xj(xo), it is easy to check 
that ([3]) and Q can be rewritten as 



hfj,{xo) — liminf 



log/i(3Aj(a;o)) 



j->+oo log 2 3 
3. A LOCAL MULTIFRACTAL FORMALISM FOR A DYADIC FAMILY 

3.1. Mult ifr act al analysis on a domain VL. 

Definition 6. Let fi be a non-empty open subset ofM. d . A collection of nonnegative 
quantities £ = (e\) indexed by the set of dyadic cubes A C fi is called a dyadic 
function on Q. 

The choice of the dyadic setting may seem arbitrary; however, it is justified by 
two reasons: 

• It is the natural choice when dealing with orthonormal wavelet bases (though 
wavelets could be defined using other division rules, in practice the dyadic 
one is the standard choice), and also the measure setting. 

• When analyzing experimental data through regressions on log-log plots, for 
a given resolution, the dyadic splitting yields the largest number of scales 
available in order to perform the regression. 

Definition 7. The pointwise exponents associated with a dyadic function £ on f2 
are the function h(x) and h(x) : Q — ¥ K defined for x £ as follows: 

• The lower exponent of £ is 

h £ (x) = lim inf l0geAj(:r) (15) 
k ; 3^+00 log 2-J v ' 

• The upper exponent of £ is 

MxHhmsup 1 " 6 ^ . (16) 
log 2 3 

By convention one sets hg(x) = hs(x) = +00 if x ^Supp(£). 

We saw in the introduction the first example of scaling function which has been 
used. We now define them in the abstract setting supplied by dyadic functions. We 
denote by A^ 2 the subset of Aj composed of the dyadic cubes contained in Q. 

Definition 8. Let fl be a nonempty bounded open subset ofM. d .The structure func- 
tion of a dyadic function £ on is defined by 

VpeM, s 3 (n 7 p) - ]T {e x y. (17) 

AeA 3 fi 

The scaling function of £ on fl is defined by 

V P eI, T?( P ) = \hvM logS ff\ (18) 

j— >+oo log Z J 
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If fl is not bounded, one defines the scaling function as follows: 

if fl n = nn 5(0, ra), Vp e R, T%(p) = Km r^{p). (19) 



> DC 



Note that the limit exists because the sequence is decreasing. From now on, we will 
assume that the set £1 is bounded, so that, at each scale j, a finite number only 
of dyadic cubes A satisfy A C !!. The corresponding results when f2 is unbounded 



follow easily from ( 19 ) 



Apart from the scaling function, an additional "global" parameter plays an im- 
portant role for classification in many applications; and, for multifractal analysis, 
checking its positivity is a prerequisite in the wavelet setting (see [2] and references 
therein): The uniform regularity exponent of £ is defined by 

log( sup e\) 

M?=liminf - AgAj . . (20) 

4 j->+oo log 2~3 X ' 

The scaling function t £ is concave (as a liminf of concave functions) taking values 
in R. The following regularity assumption is often met in practice, and implies that 
r]£ is finite for any value of p. 

Definition 9. A dyadic function £ is regular in CI if 
3d, C 2 > 0, 3A, B e R VA C Q : e A + => C X T B ^ < e x < C 2 2~ Aj . (21) 
The existence of A is equivalent to the condition hQ > —oo. More precisely, 



h £ =sup{A: the right hand side of (21) holds} 



In the measure case and in the Holder exponent case, one can pick A = 0. In the 
Holder case, the uniform regularity assumption means that A > 0. When the e\ 
are wavelet leaders, the assumption on the lower bound implies that the function 
/ considered has no C°° components. 

Since the scaling function is concave, there is no loss of information in rather 
considering its Legendre transform, defined by 

C$(H):=mHp-T$(p)). (22) 

The function £^(H) is called the Legendre spectrum of £ . 

Though it is mathematically equivalent to consider £^(H) or Tg(p), one often 
prefers to work with the Legendre spectrum, because of its interpretation in terms 
of regularity exponents supplied by the multifractal formalism. 

Definition 10. Let £ be a dyadic function on £1, and define, for H G [— oo,+oo], 
the level set associated with £ 

Eg(H) = {xeCl: h £ {x) = H}. 

4? : H e R^ dim E%(H). 

Let us now show how a heuristic relationship can be drawn between the multi- 



fractal and the Legendre spectra. The definition of the scaling function ( 18 1 roughly 
means that, for j large, Sj(Q,p) ~ 2~ Te ^ 3 . Let us estimate the contribution to 
Sj(tl,p) of the dyadic cubes A that cover the points of E £ (H). By definition of 
E £ (H), they satisfy e A - 2~ Hi ; by definition of d%(H), since we use cubes of the 
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same width 2 3 to cover E, we need about 2 d £ such cubes; therefore the cor- 
responding contribution is - 2 d " { - H ^2- Hp: > = 2-J'( ff P- d 2(ff)). When j -> +oo, the 
smallest exponent brings an exponentially dominant contribution, so that 

T2(p)=w£(Hp-d%(H)). (23) 

ri 

This formula can be interpreted as stating that the scaling function is the Legendre 



transform of the spectrum. Assuming that d^(H) is concave, it can be recovered 
by an inverse Legendre transform, leading to 

d%(H) = w£(Hp- 4\p)). (24) 

When this equality holds, the dyadic function £ satisfies the multifractal formalism, 
which therefore amounts to state that the Legendre spectrum coincides with the 
multifractal spectrum. 

Note that the derivation we sketched is not a mathematical proof, and the de- 
termination of the range of validity of ( p4| (and of its variants) is one of the main 
mathematical problems concerning multifractal analysis. The only results which 
hold in all generality are upper bounds of dimensions of singularities. 

Proposition 4. [201 [Ml [13] Let £ be a dyadic function on ft. Then 

d%{H) < C^(H). (25) 

An important consequence of this corollary is supplied by the only case where the 
knowledge of the scaling function is sufficient to deduce the multifractal spectrum, 
and even the pointwise exponent hg everywhere. 

Corollary 2. Let £ be a dyadic function. If its scaling function Tg satisfies 

3a > such that Vpel, Tg (p) = t £ (0) + ap, (26) 
then the multifractal formalism is satisfied on O, and the lower exponent of £ is 

\fx e Supp £, hs(x) = a. 



Proof, (of Corollary [2]) Assume that (26 1 is true. Then Lz(H) = — oo except for 
H = a; Corollary [2] implies in this case that d^^H) < — oo for H ^ a. Therefore 
only one Holder exponent is present, so that Vx, h(x) = a; it follows that d^?(a) = 1, 
and the multifractal formalism therefore holds. □ 

This corollary has direct implications in modeling: Indeed, several experimental 
signals have a linear scaling function. In such situations, multifractal analysis yields 
that the data have a constant pointwise exponent; therefore it supplies a non- 
parametric method which allows to conclude that modeling by, say, a fractional 
Brownian motion, is appropriate (and the slope of the scaling function supplies the 
index of the FBM), see e.g. [2] where one example of internet traffic data is shown. 
We will also see a local version of Corollary [2] which has implications in modeling: 
Corollary [5j 

3.2. Local multifractal formalism. 

Definition 11. Let £ be a dyadic function on O. The local multifractal spectrum 
of £ is the function defined by 

WH, Vx e ft, d e (x,H) = dim{E £ (H),x) (= lim df {x/r) (H)) . (27) 
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The following result, which is a direct consequence of Proposition [2j shows that 
the local spectrum allows to recover the spectrum of all possible restrictions of £ 
on a subset u> G fl. 

Corollary 3. Let £ be a dyadic function on f2. Then for any open set lo a tt, 

MH G K, d%(H) = sup d e (x,H). (28) 

Definition 12. A dyadic family £ is said to be homogenously multifractal when 
the local multifractal spectrum d £ {x, ■) does not depend on x, i.e. 

Vx G fl, Vi? G M, d £ (x,H) = d £ 1 (H). 

A local scaling function can also be defined by making the set shrink down to 
the point xq. 

Definition 13. Let £ be a dyadic function on fl. The local scaling of £ is the 
function defined by 

WH, Vx G n, T £ (x,p) = lim t> {x ' r \p). (29) 

r— >0 



Note that the right-hand side of (29) is a decreasing function of r, and therefore 
it has a limit when r — > 0. Similarly as in the multifractal spectrum case, a 
straightforward compacity argument yields that the scaling function on any domain 
u! can be recovered from the local scaling function. 

Corollary 4. Let £ be a dyadic function on f2. Then for any open set uj c 0, 

VH G E, Tg(p) = inf T £ (x,p). (30) 

Definition 14. The scaling function of a dyadic family £ is said to be homogenous 
when the local scaling function t £ (x, •) does not depend on x. 

The upper bound supplied by Corollary [4] holds for any given ball B(x, r). Fixing 
x G SI and making r — ► 0, we obtain a following local version of this result: 

Vx G n, ViT, d £ (ar, ff) < inf - T £ (x,p)) . (31) 



We will say that the multifractal formalism holds locally at x whenever (31) is an 
equality. 

As above, this result has an important consequence: In some cases, it allows 
to determine the regularity exponent at every point, even in situations where this 
exponent is not constant. 

Corollary 5. Let £ be a dyadic function. If there exists a function a:li->R such 
that the local scaling function t £ satisfies 

Vx G fl, Vp G K, t £ (x, p) = t £ (x, 0) + a(x)p, (32) 

then the multifractal formalism is locally satisfied on fl, and the lower exponent of 
£ is 

VaieO, h £ {x)=a{x). (33) 

This result is a direct consequence of (31 1 and Corollary [2] Indeed, if (32 ) holds, 
then (31) implies that d £ (x,H) = — oo if H ^ a(x). We pick now an H ^ ct(x); 
recall that d £ (x,H) — lixOr^o d £ (H); therefore 3R > such that W < r, 
d £ r \H) = — oo. In particular, H is not the pointwise exponent at x. Since this 
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argument holds for any H ^ a(x), p3| holds, and Corollary [5] follows. 

We will see an application of Corollary[5]concerning the multifractional Brownian 
Motion in Section 5.1 Combining ( pH) with Proposition [3j yields the following 



upper bound. 

Corollary 6. Let £ be a dyadic function on fl; for any open set w C fl, 

ViJ, d%(H) < sup inf (Hp - T £ (x,p)) . (34) 

It is remarkable that, though this result is a consequence of Corollary|4j it usually 
yields a sharper bound. Indeed, assume for example that the multifractal formalism 
holds for two separated regions Wi and W2 yielding two different spectra d\(H) 



and d2(H); then (34) yields max(di(H),d2(H)) whereas the global multifractal 



formalism applied to f2 = W1UW2 only yields the concave hull of max(di(-ff), d,2(H)). 



More generally, each time (34) yields a non-concave upper bound, it will be strictly 



sharper than the result supplied by Corollary |4j 

The uniform regularity exponent also has a local form: 
Definition 15. The local exponent associated with £ is the function 

hsix) = lim hf^ x ' r \ 

r->0 £ 

Note that the most general possible local exponents are lower semi-continuous 
functions, see [4"2"| . 

It would be interesting to obtain a similar characterization for the functions 
(x,H) — ¥ ds(x,H) and (x,p) — ¥ Tg(x,p) (considered as as functions of two vari- 
ables) and determine their most general form. 

3.3. An example from ergodic theory. Let £1 = (0,1). Consider a 1-periodic 
functions (f> : K — > K, as well as two continuous functions 7 : [0, 1] —> (0, 00) and 
9 : [0, 1] — s- E. Let T:i6l4 2x. For x e R and j e N denote by Sj^x) the j th 
Birkhoff sum of cj) at x, i.e., 

i-i 

S j( f>(x) = Y,cf>(T k x). 

fc=0 

Then, for any dyadic subinterval A of of generation j, let 

eA = supe-^ S ^^-^. 

When the functions 7 and 9 are constant, the multifractal analysis of the dyadic 
family £ = (e\)\ <1 Q i reduces to that of the Birkhoff averages of j<fi + &> since 
liminf^oo ^j^-*' = H if and only if liminf^oo Sj(x)/j = Hlo e^- . This is 
a now classical problem in ergodic theory of hyperbolic dynamical systems, which 
is well expressed through the thermodynamic formalism. The function log(2)r? is 
the opposite of the pressure function of — (7^ + 9), that we denote by P 1 ^(q), i.e. 



log(2)7f(p) = P 7 , e (p) = lim ilog V (supe- 7S ^ (K) - J ' e ) P Op € 

j— >oo 7 £ — ' 7 f i 
J xeAf e 

= P(- 1P )-9p, 
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where P = P_i,o is the pressure function of </>; and the following result follows for 
instance from [25] . 

Theorem 1. Let H G R; i/ien E§(H) =^ i/ and onZy i/P belongs to the interval 
[(tP )'(oo), (tJ?)'^ 00 )] a « d </l * s case if(H) = inf { Hp - rj? (p) : P eK}. 

Continuing to assume that 7 and 9 are constant, and using the fact that £ 
possesses the same almost multiplicative properties as weak Gibbs measures (see 
[101 I3T] for the multifractal analysis of these objects), i.e. some self-similarity 
property, it is easily seen that we also have Tg = rp and d£ ~ dp for all open 
subsets of il. 

Now suppose that 7 or 9 is not constant. Such a situation should be seen locally 
as a small perturbation of the case where these functions are constant, and it is 
indeed rather easy using the continuity of 7 and 9 to get the following fact. 

Proposition 5. VieH, Vg G BL 

T£[X ' P) - log(2) " ■ (35) 

Suppose also that 4> is not cohomologous to a constant, i.e. the pressure func- 
tion P of </> is not affine, which is also equivalent to saying that the interval 
/ = [P'(— 00), P'(oo)] of possible values for liminfj_ i . 00 Sj(y)/j, is non trivial. 

For all PGM, define 

Plog(2)-%) 



£h : y G (0, 1) M- 



7(2/) 



Notice that limmfj_>oo J ^-f = H if and only if liminfj-^oo Sj(y)/j = h and 
P = (7(y)/i + %))/log(2),i.e. h = £ H (y). 
Now fix a; G (0, 1). For r > we thus have 

E* {x ' r) {H) = {y G P(x,r) : liminf = frfo)}, (36) 



and due to Theorem 2.3 in [18 , for all H > 0, 

dimPf {x ' r) (H) > sup{inf {P(p) - pa : p G R} : a G rg(^|s(x,r)) n mt(P)}. 

Fix P G (t£(x,oo),t£(x,-oo)) = ( 7 (a;)P'(-°°) + 9(x),j(x)P / (00) + 6{x)). By 
construction, 

tH\ B ( x ,r)(?) = (H\og(2) - e(x))ft(x) G rg(£tf| S(:C)r) ) fl int(/). 

Thus, due to (36 1, 

dimP| (a:,r) (P) > inf |P(p) -p(Plog(2) - 0(z))/7(a:)) :p G M}, 

which, due to (35), is exactly inf{Pp — Tg{x,p) : p G R}. Since this estimate holds 
for all r > 0, 

d e (x,H) > inf {Hp -T £ (x,p): peR}, 

hence, by (31), it follows that 

d £ (x, H) = inf {Hp - T £ (x,p) :peM.}. 

For the case where H G {t^(x, 00), t' £ {x, — 00)}, it is difficult to conclude in full 
generality. We thus have proved the following result. 
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Theorem 2. Suppose <j) is not cohomologous to a constant. Fix x £ f2 and 
ff £l. IfH<£ [t^(x,oo),t^(x,-oo)} = [ 7 (a;)P'(-oo) + 9(x), 1 {x)P'(oo) + 0{x)} 
then Eg^ a ' (H) = for r small enough, and if H G (t' £ (x, oo), t^(x, — oo)) i/ien 
<fe(as, ff) = inf{ffp - t £ (x, P ) : p € K}. 

Let us mention that if the union of the sets of discontinuity points of 7 and 9 
has Hausdorff dimension 0, then the study achieved in [18 shows that the previous 
result holds at any point x which is a point of continuity of both 7 and £. Also, 
when (f> and 9 are positive, the family £ can be used to build wavelet series whose 
local multifractal structure is the same as that of £ . 



4. Measures with varying local spectrum 

4.1. General considerations. Let [i be a positive Borel measure supported by 
[0, l] d . recall that one derives from /j, the dyadic family £ fJi = {e\ := /i(3A)}AeA- 

It is obvious that the definition Q of the local dimension /i M (xo) is equivalent 
to (151 with the dyadic family £^. Similarly, the classical formalism for measures 



on [0, l] d is the same as the one described in the previous section for the family 
on ft = [0, l] d . Hence one can define a local multifractal spectrum for measures by 
Definition [271 

In the measure setting, the following result shows that the mass distribution 
principle has a local version. 

Proposition 6. Let [i be a Radon measure, A C K d and x G An supp (//). Then 

dim(x, A) > h^(x). 



Proof. It follows from (54 1 applied on AU B(x, r), remarking that the hypothesis 
x 6 supp (/x) implies that fi(A U B(x, r)) > and then letting r — > 0. □ 



We introduced the local multifractal spectrum to study non-homogeneous multi- 
fractal measures. It is interesting to recall the result of [21], where it is proved that 
homogeneous multifractal measures and non-homogeneous multifractal measures 
do not exhibit the same multifractal properties. 

Theorem 3. Consider a non-atomic homogeneous multifractal measure supported 
on [0,1]. Then the intersection of the support of the (homogeneous) multifractal 
spectrum of d M with the interval [0, 1] is necessarily an interval of the form (a, 1] 
or [a, 1], where < a < 1. 

This is absolutely not the case for non-homogenouely multifractal measures: 
consider for instance two uniform Cantor sets Cq and C\ of dimension 1/2 and 
1/4 on the intervals [0, 1/2) and [1/2, 1]. Then the barycenter of the two uniform 
measures naturally associated with Co and C satisfies 

( 1/4 if h = 1/4, 
d^h) = I 1/2 if h =1/2, 
[ —00 else. 

Hence the local spectrum is the natural tool to study non-homogeneous multifractal 
measures. 
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4.2. A natural example where the notion of local spectrum is relevant. 

The Bernoulli (binomial) measure is perhaps the most natural and simple multi- 
fractal object, and it is now folklore that is is homogeneously multifractal. We 
make a very natural modification in its construction, which will break homogene- 
ity by making the Bernoulli parameter p depend on the interval which is split in 
the construction. Doing this, we obtain a "localized" Bernoulli measure whose local 
spectrum depends on x. This example is closely related with the example developed 
in Section [331 



Let p = [0, 1] i-> (0, 1/2) be a continuous mapping. For n > 1, (ei, e 2 , e„) € 
{0, 1}, we denote the dyadic number fc eie2 ... era = Ylj—i an d the dyadic interval 

Ie 1 e 2 ...e n = [k ei e 2 ...e n , k ei e 2 ...e n +2""), where n > 1, (ei , e 2 , . .. , €„) <E {0,1}, and we 
will use the natural tree structure of these intervals using the words (ei e 2 . • ■£«)■ 

Consider the sequence of measures {nn)n>i built as follows: 

• fix is uniformly distributed on I and I\, and Hi(Io) — p(2~ r ) and /Zi(/i) = 

• H2 is uniformly distributed on the dyadic intervals I ei e 2 °f second genera- 
tion, and 

Ma(-feio) = MlC-feJ -P(fceil) and ^{hil) = Ml(-^ei) • (1 ~ p(fc ei l)). 

• ... 

• fi n is uniformly distributed on the dyadic intervals I ei e 2 ...e„ of generation 
n, and 

^n{h 1 e 2 ...e n -iOJ = Mn-l Ujeie 2 -e»-i ) ' R^Ieie 2 ...e n _il) 
and Mn(-^e 1 e 2 ...e„_ 1 l) = Mn-l(^i) ' (1 ~ K*eie a ...e„_ix))- 

Observe that by construction, for every n, for every p > n and every dyadic 
interval / of generation n, one has fJ- p (I) = /i„(J). 

Definition 16. The sequence of measures (/U n ) n >i converges weakly to a measure 
fj, that we call the "localized" Bernoulli measure associated with the map p. 

Obviously, if p is constant, one recovers the usual Bernoulli measure with pa- 
rameter p. 

We indicate the sketch of the proof to obtain the local multifractal properties of 
/i. We do not use exactly the exponent defined by Q, for simplicity we work 
with the dyadic local exponent defined by 



h d (x) = lim inf 



log H(lj{x) 



^ ' i^+oo log 2~i 

where (as usual) Ij(x) stands for the unique dyadic interval of generation j con- 
taining x. What we are going to prove also holds for the exponent /i„, but would 
require long technical developments. In particular, we would need an extension 
of Corollary 2 of [TH] on localized multifractal analysis of Gibbs measures. This 
exponent can also be encompassed in the frame of Section [3] by using the dyadic 
family £ — {^(A)}aga, thus all the "local" notions we introduced hold for this 
exponent. 

Theorem 4. For every x € [0, 1], the local spectrum associated with the exponent 
hfc of n at x is that of a Bernoulli measure of a parameter p(x) , i.e. 

Vff>0, d ll (x,H)=dn pW (H). 
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For every x *E [0,1], we consider its dyadic decomposition x — £162 — ^n---? €E 
{0, 1}. Let N 0>n {x) = #{1 < k < n : e k = 0} and N 1>n (x) = #{1 < k < n : e k = 
1} (= n — No yTl (x)). We consider the asymptotic frequencies of 0's and l's in the 
dyadic decomposition of x defined as 

N (x) = limsup -N ^ n (x). 

n— f+OO Tl 

Proposition 7. For every x G [0, 1], we have 

h%x) = -N {x)\og 2 p{x) - (1 - N (x))log 2 (l-p(x)). 

Essentially, the localized binomial measure looks locally around x like the bino- 
mial measure of parameter p(x). 

Proof. Let us fix q £ (0,1/2), and consider the classical Bernoulli measure fi q of 
parameter q on the whole interval [0, 1]. It is classical that the Holder exponent of 
\x q at every point x is 

KM) = - N °W lo S2 1 ~ (1 - W) log 2 (l - <?)• (37) 

Inspired by this formula, a Caesaro argument gives the proposition. Indeed, by 
construction, the value of the /i-mass of the interval /„ (a;) is given by 



fl(I n (x)) = JJp(fc eie2 ...e»_ 1 l)*, 



i=l 



where 



P{ke 1 e 



A* = )P( k ^2-ei-ii) ife 4 = 

|l -p(fc ei e 2 ...e j _ 1 l) if 6, = 1 

Hence, 

fi(I n (x)) — 2^'*= 1:e i= los 2P( fe n=2- -^-ii)+Er = i :ei= o ^(i-pCfen^- ^i-ii)) 

Since p(A; 6l£2 ... 64 _ 1 i) tends to p(x) when i tends to infinity, and since N (x) is 
the asymptotic frequency of zeros in the dyadic expansion of x, one sees that 

1 " 

- zl log 2 p(fc ei e 2 ... e^l) > n -, +00 N (x)p(x). 

n * — ' 

j=l: ei =0 

Similarly, since p[x) < 1/2, 
1 - 

- >J l0g 2 (l -K fc ei£2...e 4 -il)) ►n-^+oo (1 ~ N {x)){l - p(x)). 

n * — ' 

i=l:e<=l 

Let a — —N (x) \og 2 p{x) — (1 — N (x)) log 2 (l — p(x)). The latter proves that, 
given e > 0, there exists an integer N such that n > N implies that 

This yields the result. □ 

Consider an interval J c [0,1], and the multifractal spectrum d^(H,J) = 
dim {a; € J : h d ^(x) — H}. The value of this spectrum is a consequence of the 
following theorem of Barral and Qu in [T5] (who proved this result for any Gibbs 
measure /it). 
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Theorem 5. Fix q G (0, 1/2), and consider the Bernoulli measure with parameter 
q. Let us denote by R q the support of the (homogeneous) multifractal spectrum of 
fi q . Let h : [0, 1] —> R q be a continuous function. Then for every interval J C [0, 1], 
one has 

dim {a; G J : /ijj (x) — h(x)} — supjd^ (h(x)) : x G J}. 
We now prove Theorem [4] 

Fix H > 0, and also g€ (0, 1/2). If for some x one has 

-N (x) log 2 p(x) - (1 - iVo(aO) log 2 (l - p(a;)) = h, 

then there exists a real number h q (H,x) such that 

-iVo(x) log 2 q - (1 - N (x)) log 2 (l - q) = h q (H, x). 

Since both p(x) and q are strictly less than 1/2, a simple argument entails that the 
map h q {H, x) is continuous with respect to {H, x). 

Now fix xq G [0, 1] and consider the Bernoulli measure with parameter q = p(xq). 
Consider the interval I = B(xo,r). One has 

{x e I : K{x) =H} = {x e I :-N Q (x)log 2 p(x) - (1 - N (x)) log 2 (l - p(x)) = H) 
= {xel: -N Q (x) log 2 q - (1 - iV (x)) log 2 (l - q) = fe,(ff, a;)}. 
But this last set has its Hausdorff exactly given by Theorem [5j hence 
dim{x € J : = if} = sup{(i M (h q (H, x)) : x G J}. 

When r goes to zero, p(x) tends uniformly to q = p(xo)- Hence h q (H,x) tends to 
H. In particular, the mapping being continuous (real analytic in fact), when r 
goes to zero one finds that 

d fi (x,H) = d llq (H) = d flpU) (H), 

This result can immediately be applied to the case where the mapping x t— > p(x) 
is continuous by part (instead of simply continuous), and can certainly be adapted 
when p is cadlag. It would be worth investigating the case where p enjoys less 
regularity properties. 

Remark 1. Many examples of Cantor set with varying local Hausdorff dimensions 
have been constructed (51151). here the key point is that we perform the (global and 
local) multifractal analysis of measures sitting on these "inhomogeneous" Cantor 
sets. 



5. Local spectrum of stochastic processes 

Suppose now that / is a nowhere differentiable function defined on [0, l] d ; one 
can associate with / the dyadic family £f = {0sc/(3A)}asA, where the oscillation 
of / over a set w C f! is 

OsCfiuS) = sup{/(x) : x G uj} — inf{/(x) : x G a;}. 

Then, it is obvious that the pointwise Holder exponent (Ifj]) of / at x is the same 
as the one defined by (15) with the dyadic family £f. Hence, the previous de- 
velopments performed in the abstract setting of dyadic functions family holds for 
non-diffcrcntiablc functions. 
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We start by giving a simple general probabilistic setting which naturally leads 
to a weak, probabilistic form of homogeneity. Let X be a random field on R d ; X 
has stationary increments if Vs £ K d , the two processes 

x h-> Y s {x) := X(s + x) — X(s) and x h-> X(x) 

share the same law. Indeed, this equality in law implies the equality in law of the 
linear forms applied to the two processes Y s and X, hence of iterated differences 
and wavelet coefficients. It follows that local suprema of iterated differences and 
of wavelet coefficients computed on dyadic cubes also share the same laws, and 
Proposition [3] implies that, if X has locally bounded sample paths, then the Holder 
exponent has a stationary law. Therefore, the Holder spectra on dyadic intervals 
of the same width also share the same law almost surely. As a result, the Holder 
spectra on all dyadic intervals share the same law. This leads to the following 
result. 

Proposition 8. Let X be a random field on R d with stationary increments. If X 
has locally bounded sample paths, then 

Vs a.s. V-ff d x (s,H) =d x (0,H). 

5.1. Local analysis of the multifractional Brownian motion. Let H denote a 
function defined on K d with values in a fixed compact subinterval [a, b] of (0, 1). We 
assume that H satisfies locally a uniform Holder condition of order f3 £ (b, 1), that 
is, H £ CP(Q) for every open subset f2 of K d . Now, recall that the multifractional 
Brownian motion (MBM) with functional parameter H has been introduced in (T9l 
149] as the continuous and nowhere diffcrentiable Gaussian random field Bh = 
{B H {x), x £ R d } that can be represented as the following stochastic integral 

Bh{x) = L **®> 

where x-f denotes the standard inner product, |f I2 is the usual Euclidean norm, and 
dW stands for the "Fourier transform" of the real-valued white noise dW , meaning 
that for any square-integrable function /, one has 

f(t;)dw(z)= ( f{x)dw{x). 



In particular, the MBM reduces to a fractional Brownian motion when the function 
H is chosen to be constant. The pointwise regularity of the MBM is well known; 
as a matter of fact, it has been shown in [7j that 

a.s. y X £R d h BH (x)=H{x). (38) 

Thus, the Holder exponent of the MBM is completely prescribed by the function 
H . Our purpose is now to give an illustration to Corollary [5] above by showing 
that the multifractal formalism is locally satisfied by almost every sample path of 
the MBM. To be specific, we shall establish in the remainder of this section the 
following result which, with the help of Corollary [5] enables one to recover ( 38 ) . 



Proposition 9. Let Eh denote the dyadic function that is obtained by considering 
the wavelet leaders of the multifractional Brownian motion Bh, and assume that the 
wavelets belong to the Schwartz class. Then, the local scaling function Tg H satisfies 

a.s. Vx£R d VpeM t £h (x, P ) = H(x)p- d. 
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In order to establish Proposition |9j we shall work with a Lemarie-Meyer wavelet 
basis of L 2 (R d ) formed by the functions 2 d ^ 2 ip^ (2^ x — k), see [H], and more 
generally with the biorthogonal systems generated by the fractional integrals of the 
basis functions namely, the functions i/jW> /i defined by 



m h 2 +d/2 ' 



It will also be convenient to consider the Gaussian field Y = {Y(x,h), (x,h) g 
R d x (0, 1)} given by 

Note, in particular, that Bh(x) = Y(x, H(x)) for all x € R d , and that the random 
field {Y(x, h), x £ M. d } is merely a fractional Brownian motion with Hurst param- 
eter h. By expanding its kernel in the orthonormal basis of L 2 (M. d ) formed by the 
Fourier transforms of the functions 2*'/ 2l 0W(2^ar— k), and by virtue of the isometry 
property, the stochastic integral defining Y(x,h) may be rewritten in the form 

Y( x , /o = EE E 4,k 2 ' hj (V^z - k) - ^ h {-k)) , 

where the k form a collection of independent standard Gaussian random variables. 
It is possible to show that the above series converges uniformly on any compact 
subset of R d x (0, 1), see [BJ. Moreover, the above decomposition yields the following 
natural wavelet expansion of the field Bh : 

5j*(a0 = EE E Ak 2 ~ H{x)j (^ {l) ' H{x) {Vx -k)- ^' H(x \-k)^j . (39) 

Furthermore, it is shown in [BJ that the low-frequency component of Y, that is, 

-l 

E E E s U 2 ~ hj (V^^z - *) - v w ' ft (-*o) , 

i j = -oo fegZ d 

is almost surely a C°° function in the two variables x and h. Hence, the low- 
frequency component of the MBM, which is obtained by summing only over the 



negative values of j in (39 1, is in C"(Cl) for any open subset f2 of M. d , just as the 
functional parameter H . As f3 is larger than all the values taken by the function 
H , it follows that the pointwise regularity of the MBM is merely given by that of 
its high-frequency component, that is, 

oo 

b h (x) = E 4fc 2 ~ H( * (^ t)M{x H^x -k)- v (i) ' H(x) (-fc)) . 

As a consequence, we may consider in what follows the high-frequency component 
Bh instead of the whole field In addition, in view of the regularity of H, it 
follows from standard results on Calderon-Zygmund operators (see [IS]) and robust- 
ness properties of the local scaling functions, T£ H coincides with the local scaling 
function of the dyadic family En which is obtained by considering the wavelet 
leaders associated with the wavelet coefficients 

c j,k - £ j.k L 
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(Recall that in [34| , it is proved that the scaling function is "robust" , i.e. does not 
depend on the smooth enough wavelet basis chosen; furthermore, the arguments of 
the proof clearly are local, so that the local scaling function also is robust.) 



Letting A denote the cube corresponding to the indices i, j and k as in Section 2.3 
these coefficients may naturally be rewritten in the form 

where s\ is the standard Gaussian random variable fc , x\ is the basis point k2~ : > of 
the cube A and (A) is its scale j. Recall that the wavelet leaders d\ are then defined 
in terms of the wavelet coefficients through (13). Finally, for the sake of simplicity 
and without loss of generality, we shall study the local scaling function Tg H only on 
the open set (0, l) d , so that we only have to consider the dyadic subcubes of [0, l) d . 

Let us now establish a crucial lemma concerning the behavior on the subcubes 
of [0, f ) d of the new dyadic family £h- 

Lemma 1. With probability one, for any dyadic cube A C [0, l) d with scale (A) 
large enough, 

(A)3H(xa) - ax - hi- 

proof. We begin by the proving the lower bound. For any proper dyadic subcube 
A of [0, l) d with scale (A) = j, we have 

¥(d x < (A) _3H(:EA) 2~- H(;rA)<A) ) = Y\ F (I £ A'| < (X)~ 3H ^ Xx) 2 H< - x ^'^ x " > ~ H ^ W ). 

A'C3A 

Let l(j) = j + [~(2/<i)log 2 f\, where [•] denotes the ceiling function and log 2 the 
base two logarithm. Considering in the above product only the subcubes A' C 3A 
with scale (A') equal to l(j), and using the elementary fact that the modulus of a 
standard Gaussian random variable is bounded above by t with probability at most 
t, we deduce that 

p(d A < ( X y3H(x x ) 2 -H(x x )(\)-j < -Q ^-3H(x x ) 2 H(x x ,)(\')-H(x x ){X) _ 

A'C3A 

Moreover, the function H satisfies locally a uniform Holder condition of order /3, 
so there exists a real C > that does not depend on A such that 

VA'c3A \H(x y ) - H(x x )\ < C2- pi . (40) 

Combined with the observation that there are at least j 2 subcubes A' C 3A such 
that (A') = l(j), this implies that 

P(^A < (A)- 3 ^^--^^^) < ^j-3H(x x ) 2 H(x x )(lU)-j)+Cl(f)2-^y 

Given that the function H is valued in the interval [a, b] , we infer that 

P(^A < (A)- 3 ^ 2 ^--^^) < ^j(2/d-3)a 2 6+Ci(j)2-«y 

■ 2 

The right-hand side is clearly bounded above by e~ J when j is larger than some 
integer jo, so that 

J2 ^ (A)- 3 ^)2- H ^)< A >) < 2 * e_ ' 2 < °°. 

AC[0,l) d j>jo 

W>j 
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and we deduce the required lower bound from the Borel-Cantelli lemma. 

In order to establish the upper bound, let us begin by observing that with proba- 
bility one, for any dyadic cube A C [0, l) d with scale (A) = j large enough, \e\\ < j. 
This follows again from the Borel-Cantelli lemma, together with the fact that 

P(| £A |>j) = 2(i-<i>(j))<^— 

J 

which itself follows from standard estimates on the asymptotic behavior of the 
cumulative distribution function $ of the standard Gaussian distribution. Now, 



along with (401, this implies that for (A) = j large enough, 



dx < sup (A')2-W**>- ca ~ M )< v > = j 2 -W**)-<n-")i < 2j2~ H( - x ^ j , 

A'C3A 

and the required upper bound follows. □ 

We may now finish the proof of Proposition [9] To this end, let x e (0, l) d and 
r > such that Q — B(x,r) C (0, l) d . Then, owing to Lemma [l] the structure 



function of the dyadic function Eh on Q,, which is defined by (17), satisfies 



' <Si(0,p)< J2 (2&- H{xx) i) (41) 



I j3H(x x 

AeA" V J AeA 



for j large enough and p > 0. Given that H satisfies locally a uniform Holder 
condition of order f3, there exists a real C > that depends on neither x nor r such 
that \H(x\) — H(x)\ < CrP for all dyadic cubes A C f2. In addition, the cardinality 
of A^ 2 is comparable with r d 2 dj . Thus, there is a constant C > such that 



r d 2 dj I 2 -(H(x)+Cr l3 )j 



C> \ j3(H(x)+Crf>) 



< Sj(tl,p) < C'r d 2 d J (2 3 2-( H ^- Cr ^y . 



It follows that the scaling function of £ h on fl satisfies 

(H(x) - Cr )p -d<rp (p) < (H(x) + Cr p )p - d. 

Letting r go to zero, we may finally conclude that T£ H (x,p) = H(x)p — d for all 
p > and x £ (0, l) d . The same approach still holds for the negative values of 



p except that the inequalities have to be reversed in ( 41 ) and in the subsequent 
estimates as well. Proposition [9] follows. 

5.2. A Markov process with a varying local multifractal spectrum. In this 
section we reinterpret the results of [3] in terms of local spectrum. A quite gen- 
eral class of one-dimensional Markov processes consists of stochastic differential 
equations (S.D.E.) with jumps. Recall that such a process is the sum of a Brow- 
nian motion and a pure jump process. We will assume in the following that the 
process has no Brownian part; indeed, since Brownian motion is mono-Holder, its 
consequence on the spectrum is straightforward to handle: it eliminates Holder ex- 
ponents larger than 1/2 and, eventually adds a point at (1/2, 1). Thus the Markov 
processes that will be studied are jumping S.D.E. without Brownian and drift part, 
starting e.g. from 0, and with jump measure v(y, du) (meaning that when located 
at y, the process jumps to y + u at rate v(y,du)). Again, since this is a "toy" 
model, we will make additional simplifying assumptions: Namely that the process 
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is increasing (that is, u(y, (— oo, 0)) = for all y 6 M). Classically, a necessary 
condition for the process to be well-defined is that J Q uv{y,du) < oo. 

If v is chosen so that the index $ v < v \ is constant with respect to y, then one 
expects that the local multifractal spectrum d,M(t,h) of the process M — (M t ) t > 
will be deterministic and independent of t. Hence, the index of the jump measure 
will depend on the value y of the process. The most natural example of such a 
situation consists in choosing 

v 7 (y,du) := 7(y)u~ 1 ~ 7(y) l [0!l] (u)cZM, 

for some function 7 : Mh> (0, 1). The lower exponent of this family of measures is 

> 0, P Vy ( y ,.) = i{y)- 

In |9_, the following assumption is made 



(H) 



{There exists e > such that 7 : [0, 00) 1 — > [e, 1 — e] 
is a Lipschitz-continuous strictly increasing function. 



It is relatively clear that the assumptions can be relaxed, and that many classes 
of Markov processes could be further studied. An interesting subject to investigate 
is the range of functions 7 that could be used in the construction. For a process, 
M = (M t ) t >o, one sets AM t =M t - M t -, where M t _ = lim M s 



¥t, S<t 



Proposition 10. [9] Assume that {%) holds. There exists a strong Markov process 
M = (M t ) t > starting from 0, increasing and cddldg (i.e. right- continuous, with 
a left limit), and with generator C defined for all y £ [0, 00) and for any function 
<f> : [0,oo) i—)- K Lipschitz-continuous by 

C<f>{y)= ( [<f>(y + u)-<f>(y)}^( y) du). (42) 



./() 

Almost surely, this process is continuous except on a countable number of jump 
times. Denote by J the set of its jump times, that is J = {t > : AM(t) 0}. 
Finally, J is dense in [0,oo). 

This representation of M is useful for its local regularity analysis. 

The following theorem of [S] summarizes the multifractal features of M. 

Theorem 6. Assume {W) and consider the process M constructed in Proposition 
\1C\ Then, the following properties hold almost surely. 

(i) For every t € (0, oo)\j7 ; the local spectrum of M at t is given by 

' [-00 if h>lMM t ), v ' 

while for t G J , 

'h"y(M t ) if < h < l/ 7 (M t ), 
d M (t, h) = {h- 7 (M t _) if h g [lh(M t ), l/7(M t _)], (44) 
-00 if h > l/7(Af t _). 
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(ii) The spectrum of M on any interval I — (a, b) c (0, +00) is 

yh>0, d M {h) = sup j/i • 7(M t ) : te I, h- j(M t ) < 1 j (45) 

= sup I h ■ 7(M S _) : sejni, h- j(M s -) < l}. (46) 



In (45) and (46), we adopt the convention that sup( 





Figure 1. Two sample paths of the stochastic process M built 
using the function 7(2/) := min(l/2 + y/4, 0.9). On the right hand- 
side are plotted the theoretical spectra Dm{[0, 3], .). 



As can be seen from the definition of the local multifractal spectrum, in order 
to prove Theorem [6j it is enough to show ( 46 1 . Indeed, ( 43 ) simply follows from 
considering the limit of (46) when the interval I is the centered ball B(t,r) and 
letting r tend to zero. 



Formula (46) is better understood when plotted: for every s G J fl J, plot a 
segment whose endpoints are (0,0) and (1/7(M S _), 1) (open on the right), and 
take the suprcmum to get Dm(I,.). Sample paths of the process M and their 
associated spectra are given in Figure 1. 

The formulae giving the local and global spectra are based on the computation 
of the pointwise Holder exponents at all times t. The value of the pointwisc Holder 
exponent of M at t depends on two parameters: the value of the process M in 
the neighborhood of t, and the approximation rate of t by the set of jumps J . In 
particular, the following properties holds a.s., 

for every t > 0, h M (t) < l/l(M t ), 
for Lebesgue-almost every t, ^jw(i) = ^/l(M t ), 

for every k € (0, 1), dim ff {< > : h M (t) = n/~f(M t )} = n. 
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The relevance of the local spectrum in this context is thus obvious: depending on 
the local value of M, the pointwise Holder exponents change, and so is the (local) 
multifractal spectrum. 

It is worth emphasizing that, as expected from the construction of the process 



M, the local spectrum (43) at any point t > essentially coincides with that of 
a stable Levy subordinator of index r y(M t ). This local comparison is strengthened 
by the following theorem, which proves the existence of tangent processes for M 
(which are Levy stable subordinators). 

Proposition 11. We denote by T t := a{{N(A),A e B{[0,t] x [0,oo))». Let 
to > be fixed. Conditionally on J-" to) the family of processes 



a i/7(M* ) /te[o,i] 

converges in law, as a — > + , to a stable Levy subordinator with Levy measure 
7(.Mt )tt -1-T (* ft o)<i'ii. Here the Skorokhod space of cadlag functions on [0,1] is 
endowed with the uniform convergence topology. 

Observe that for all s € J, all h E (l/j(M s ), l/j(M s -)}, d M (h) = h ■ j(M 3 -). 
Thus the spectrum du of M on an interval / is a straight line on all segments of 
the form (l/^f(M s ), 1/7(M S _)], s E JT\I. By the way, this spectrum, when viewed 
as a map from R + to R + , is very irregular, and certainly multifractal itself. This 
is in sharp contrast with the spectra usually obtained, which are most of the time 
concave or (piecewise) real-analytic. Hence, the difference between the global and 
the local multifractal spectra is stunning: While g?m is very irregular, dM{t, ■) is a 
straight line. 

This example naturally leads to the following open problem, which would express 
that a natural compatibility holds for local multifractal analysis: Find general 
conditions under which a stochastic process X which has a tangent process at a 
point xq satisfies that the multifractal spectrum of the tangent process coincides 
with the local spectrum of X ah xq. 

6. Other regularity exponents characterized by dyadic families 

Other exponents than those already mentioned fit in the general framework given 
by Definition [7] and therefore the results supplied by multifractal analysis can be 
applied to them. We now list a few of them. 

Pointwise Holder regularity is pertinent only if applied to locally bounded func- 
tions. An extension of pointwise regularity fitted to functions that are only assumed 
to belong to L v loc is sometimes required: The corresponding notion was introduced 
by Calderon and Zygmund in 1961, see [52], in order to obtain pointwise regularity 
results for elliptic PDEs. 

Definition 17. Let p E [1, +oo) and a > —d/p. Let f E L? oc (Q), and xq E f 
belongs to T^(xq) if there exist C > and a polynomial P of degree less than a 
such that, for r small enough, 

( i / \f(x)~P(x-x )\ p dx) <Cr a . (47) 

The p-exponent of f at xq is 

h p f (x )=sup{a:feTV(x )}. 

Remarks: 
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• The normalization chosen in (47) is such that cusps \x — xo\ a (when a ^ 2N) 
have an Holder and a p-exponent which take the same value a at Xq. 

• The Holder exponent corresponds to the case p = +00. 

• We only define lower exponents here: Upper exponents could also be defined 
in his context, by considering local LP norms of iterated differences. 

• Definition [XT] is a natural substitute for pointwise Holder regularity when 
functions in L v loc are considered. In particular, the p-exponent can take 
negative values down to —d/p, and typically allows to take into account 
behaviors which are locally of the form 

l ^—rz for 7 < d/p, (48) 

\x - x p 

A pointwise regularity exponent associated with tempered distributions has been 
introduced by Y. Meyer: The weak scaling exponent (see [46] . and also [1] for a 
multifractal formalism based on this exponent). It can also be interpreted as a 
limit case of other exponents for distributions, which can be related with the Holder 
exponent; let us briefly recall how this can be done. 

Let / be a tempered distribution defined over K d . One can define fractional 
primitives of order s of / in the Fourier domain by 

/^ 7) (o = (i+ia 2 r /2 /(o- 

Since / is of finite order, for s large enough, /(- s ) locally belongs to LP (or L°°). 
It follows that one can define regularity exponents of distributions through p- 
cxponents (or Holder exponents) of a fractional primitives of large enough order. If 
/ is only defined on a domain Jl, one can still define the same exponents at xq € 17 
by using a function g e V(R d ) such that g is supported inside Q and g{x) = 1 in 
a neighborhood of x ; then fg is a tempered distribution defined on WL d and the 
exponents of (fgY~ 8 ^ at xq clearly do not depend on the choice of g. 

Let / be a tempered distribution defined on a open domain. Denote by h s Ax) 
the Holder exponent of (which is thus canonically well defined for s large 

enough. By definition, the weak scaling exponent of / at x is 

W f (x)= lim (h' f (x)-s) 

(note that the limit always exists because the quantity considered is an increasing 
function of s). We will not deal directly with this exponent because it does not 
directly fit in the framework given by Definition [7J But we will rather consider the 
following intermediate framework. 

Definition 18. Let f be a tempered distribution defined on a non-empty open 
set SI C R d . Let p > 1 and s be large enough so that /'~ s ) belongs to L p in a 
neighborhood of Xq. The fractional p- exponent of order s of f at xq is defined by 

h p f s (x ) = h p f( _ 3) (x ) 

(using the convention hy =hf). 

Note that, in practice, the standard way to perform the multifractal analysis 
of data that are not locally bounded is to deal with the exponent h°^ ,s , where s 
is chosen large enough so that € L^ c , i.e. it consists in first performing 

a fractional integration, and then a standard multifractal analysis based on the 
Holder exponent, see [2] and references therein. 
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Similarly, in the function case, if the pointwise regularity exponents are small 
enough, they can be recovered for the oscillation of /. Recall that the oscillation of 
/ of order I on a convex set A is defined through conditions on the finite differences 
of the function /, denoted by Aff f: The first order difference of / is 

(Alf)(x) = f(x + h)-f(x). 

If I > 1, the differences of order I are defined recursively by 

(Kf)(x) = (Aj^/Xz + h)- (A^fXx). 

Then 

Osc l f (A) = sup |K/)(z)|. 

x,x-\-lh£ A 

One easily checks that the Holder exponent can be derived for the oscillation on 
the cubes 3A. Let / be locally bounded on an open set fi. 
If I > hf(xo), then 

Vx G Q, h f (x ) = liminf /v " . (49) 

j— >+oo log Z J 

Recall also Proposition[3]which allows to derive numerically the Holder exponent 
by a log- log plot regression bearing on the the d\.(x ) when j — > +oo, see [3^| . 

However, in contradistinction with the measure case, a similar formula does not 
hold for the upper Holder exponent, see [23] where partial results in this direction 
and counterexamples are worked out. 

We now turn to the wavelet characterization of the p-exponent. We will assume 
that / locally belongs to L p , with slow i p -increase, i.e. satisfies 

3C,7V>0 f \f(x)\ p dx <C{1 + \R\) N . 

In the following, when dealing with the T p regularity of a function /, we will 
always assume that, if / is defined on an unbounded set 17, then it has slow L p - 
increase, and, if ^ M. d , then the wavelet basis used is compactly supported. 

Definition 19. Let f € L p oc (fl), and let %p\ be a given wavelet basis. The local 
square function of f is 

SM = { E I c v| 2 lv(a0) » 

\A'C3A / 

and the p-leaders are defined by d? x = 2 d ^ p \\ Sj.\ \\ p . 

The following result of yields a wavelet characterization of the p-exponent 
which is similar to ( 14 ) . 

Proposition 12. Let p G (1, oo) and f G L p . Assume that the wavelet basis used 
is r -smooth with r > ^(xq) + 1. Then 

logeff , s 

h> f (x )= liminf A ff . (50) 

J j->+oo log Z J 
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Recall that the "almost-diagonalization" principle for fractional integrals on 
wavelet bases states that, as regards Holder regularity, function spaces or scal- 
ing functions, one can consider that a fractional integration just acts as if it were 
diagonal on a wavelet basis, with coefficients 2~ SJ on This rule of thumb is 
justified by the fact that a fractional integration actually is the product of such a 
diagonal operator and of an invertible Calderon-Zygmund operator A such that A 
and A~ x both belong to the Lemarie algebras M? , for a 7 arbitrarily large (and 
which depends only on the smoothness of the wavelet basis) see [HI |33] for the defi- 
nition of the Lemarie algebras and for the result concerning function spaces and [34] 
and references therein for Holder regularity, function spaces or scaling functions. 

It follows from Proposition [12] , and the "almost-diagonalization" principle for 
fractional integrals on wavelet bases, that the exponent h P f S {x$) can be obtained 
as follows. 

Corollary 7. Let p £ (l,oo) and f £ L p . Let 

S h( x ) = [ E la-'cvfMao) and d% s = || S ' J>X \\ 

\A'C3A / 

Then, if the wavelet basis is r-smooth with r > h P j{xo) + s + 1, then 

logdw > 

1 ]~>+oo log 

7. A FUNCTIONAL ANALYSIS POINT OF VIEW 

7.1. Function space interpretation: Constant regularity. If p > 0, the scal- 
ing function has a function space interpretation, in terms of discrete Besov spaces 
which we now define. Recall that the elements of a dyadic family are always non- 
negative. 

Definition 20. Let set and p E R. A dyadic function £ belongs to bp'°°(f2) if 
3C, Vj, 2-*' J2 *( e A) p < C ■ T s ™. (52) 

If p — +00, a dyadic function £ belongs to b^°°(f2) if 

3C VA : e x <C-2- sj . (53) 

Note that, if p > 0, this condition (if applied to the moduli of the coefficients) 
defines a vector space. It is a Banach space if p > 1, and a quasi-Banach space if 
< p < 1; recall that, in a quasi-Banch space, the triangular inequality is replaced 
by the weaker condition : 

3C, Vx,y, || z + 2/ ||<C(|| x || + 112,11). 

Definition [20] yields a function space interpretation to the scaling function when 
p > 0. It is classical in this context to rather consider the scaling function 

Vs(p) = T £ {p) - d. 

Then, if is a bounded set, 

Vpel, r,2(p)=sup{s:£eh s p ^ ao (fl)}; 
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and, if O is unbounded, then the function space interpretation is the same, using 



the precaution supplied by (19). Additionally, 

fc£ = sup{«:£€b£»(fi)}. 
The terminology of "discrete Besov spaces" is justified by the fact that, if the 



e\ are wavelet coefficients, then |52|) and (53) are the wavelet characterization of 
the "classical" Besov spaces 13^°° (R d ) of functions (or distributions) defined on K d ; 
therefore each wavelet decomposition establishes an isomorphism between the space 
b^°°(E rf ) and the space S^°°(K d ), see [44]. Note that, when p = oo, these Besov 
spaces coincide with the Holder spaces C s (IR <i ), so that, when the (e>) are wavelet 
coefficients, then the uniform regularity exponent has the following interpretation 

h% = sup{s:£ EC s {M. d )}. 

In the measure case, H. Triebel showed that the discrete Besov conditions bearing 
on the /i(3A) can also be related with the Besov regularity of the measure /i, see 

If s < d, (m(3A)) G b;'°°(m d ) <s=> fx G B s p - d '°°(R d ). 

In the case where p = +oo, uniform regularity gives an important information 
concerning the sets A such that p(A) > 0, as a consequence of the mass distribution 
principle, see Section [2j Since this estimate precisely means that the sequence 
( e \) — (a*(3A)) belongs to b^°°(f2), it follows that, if A C £1 and if a measure /x 
satisfies p(A) > 0, then 

dim(A) > h%. (54) 

When the sequence £ is composed of wavelet leaders, or of p-leaders, the cor- 
responding function spaces are no more Besov spaces, but alternative families of 
function spaces, the Oscillation Spaces, see [331 135] . 



h% = sup{A : Q holds}. 

The following upper bounds for dimensions are classical for measures, see [20] . 
and are stated in the general setting of dyadic functions in [37] . 

Proposition 13. Let £ be a dyadic function, and let 

J% = {x G fl : h £ (x) >H}, G^{i£H: h £ (x) > H}, 

F% = {xen: h £ (x) < H)}, K% = {x G Q : h £ (x) < H)}. 

• If£e b s p '°°(Cl) with p>0, then dim(Gg) < d - sp + Hp. 

• If££ b s p '°°{n) with p>0, then dim p (F^) <d- sp + Hp. 

• If£e b*'°°(0) with p<0, then dim(JTg) < d- sp + Hp. 

• If£e %°°(r2) with p<0, then dim p ( jg) < d- sp + Hp. 

7.2. Function space interpretation: Varying regularity. Recall that the global 
scaling function has a function space interpretation in terms of Besov spaces which 
contain the dyadic function £ . Similarly, the local scaling function can be given 
two functional interpretations; one is local, and in terms of germ spaces at a point, 
and the second is global, and is in terms of function spaces with varying smooth- 
ness. We now recall these notions, starting with germ spaces in a general, abstract 
setting. 
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Definition 21. Let E be a Banach space (or a quasi-Banach space) of distributions 
satisfying D <-} E <—} D' . Let x £ Mr; a distribution f belongs to E locally at x 
if there exists tp € D such that ip(x) = 1 in a neighbourhood of x and ftp € E. We 
also say that f belongs to the germ space of E at x, denoted by E x . 

Let us draw the relationship between the local scaling function and germ spaces: 
If the (e\) are the wavelet coefficients of a function /, then 

Vp > 0, Vf (x,p) = sup [s-.fe b;Ip^] . 

Note that, in the wavelet case, these local Besov regularity indices have been inves- 
tigated by H. Triebel, see Theorem 4 of [BU] where their wavelet characterization is 
derived, (the reader should be careful that what is referred to as "pointwise regular- 
ity" in the terminology introduced by H. Triebel is called here "local regularity"). 
The uniform exponent can also be reformulated in terms of of Holder spaces: 

H f (x) = sup {s:f€C*}. 

In that case, the function 7if(x) is called the Local Holder exponent of /. Its 
properties have been investigated by S. Seuret and his collaborators, see e.g. [42] . 

Note that the definition of germ spaces can be adapted to the dyadic functions 
setting. 

Definition 22. Let E be a Banach space (or a quasi-Banach space) defined on 
dyadic functions over Q ; a dyadic function (e\) belongs to E x if there exists a 
neighbourhood w of x such that the dyadic function (ex) restricted to to belongs to 
E. 

It is clear that this definition, when restricted to the case of Besov spaces and 
wavelet coefficients coincides with Definition I2T1 

We now turn to function spaces with varying smoothness. Such spaces were 
initially introduced by Unterberger and and Bokobza in [611 162j , followed by many 
authors (see [52] for an extensive review on the subject). A general way to introduce 
such spaces is to remark that the classical Sobolev spaces H s ' p (R d ) can be defined 
by the condition 

II T(f) || p <oo, 
where T is the pseudo-differential operator defined by 

(Tf)(x) = ( e**(l + ie| 2 ) s/2 /(£R. 

This definition leads to operators with constant order s because the symbol (1 + 
|£| 2 ) s / 2 is independent of x. However, one can define more general spaces, with 
possibly varying order if replacing (1 + |£| 2 ) s / 2 by a symbol a(x, £). In particular 
the symbols (1 + |£| 2 )«0)/ 2 will lead to Sobolev spaces of varying order H a ' p where 
we can expect that, if a is a smooth enough function (say continuous), then the 
local order of smoothness at x will be a{x). This particular case, and its extensions 
in the Besov setting, has been studied by H.G. Leopold, followed by J. Schneider, 
Besov, H. Triebel, A. Almeida, P. Hasto, J. Vybfral, and several other authors, 
who gave alternative characterizations of these space in terms of finite differences 
or Littlewood-Paley decomposition. They also studied their mutual embeddings 
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(and also in the case where both the order of smoothness and the order of in- 
tegrability p vary) and their interpolation properties, see J53J [S3] fand references 
therein, and also [55] for an historical account. The reader can also consult [51 155] 
for recent extensions in particular when both the order of smoothness and the order 
of integrability p vary. We follow here the presentation of J. Schneider, since this 
author obtained Littlewood-Paley characterizations, which are clearly equivalent to 
the wavelet characterization that we now give. For the sake of simplicity, we as- 
sume form now on that the distributions considered are defined on E d and that the 
wavelet basis used belongs to the Schwartz class (the usual adaptations are stan- 
dard in the case of functions on a domain, or for wavelets with limited regularity). 
We additionally assume that the function a is uniformly continuous and satisfies 

3c, C > 0, Vx e R d , c< a{x) < C. (55) 

Then the Besov space B^ ,q (for p,q £ (0, oo]) can be characterized by the following 
wavelet condition, which is independent of the wavelet basis used. 



Proposition 14. Let a be a uniformly continuous function satisfying (55), and 
let p,q € (0, oo]. The Besov space of varying order B^ q is characterized by the 
following condition: 

Let C\ denote the wavelet coefficients of a distribution f , and let 

2"*' (cA2 a(A)j ) p I , 
XeA s J 

where a(\) denotes the average of the function a on the cube X; then f £ B p ' q if 

(aj) el". 

Note that when p = q = 2 one recovers the Sobolev space H a ' 2 defined above, 
and when a is a constant equal to s, then one recovers the standard Besov space 
Bp' 9 . Furthermore, the embeddings 



p 7 p 



yield easy to handle "almost characterizations" of Sobolev spaces of varying order. 
The following result, which follows directly from the definition of the local scaling 



function (Definition [8]) and the characterization supplied by Proposition 14 gives 
the interpretation of the local scaling function in terms Besov spaces of varying 
order. 

Proposition 15. Let f be a distribution defined on M. d . Then, for p > 0, the local 
wavelet scaling function of f can be recovered by 

Vp > 0, n f (p, x)=p- sup{a : / € B a p °°}. 
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